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ABSTRACT 



We, first, analytically work out the long-term, i.e. averaged over one orbital 
revolution, perturbations on the orbit of a test particle moving in a local Fermi 
frame induced therein by the cosmological tidal effects of the inhomogeneous 
Lemaitre-Tolman-Bondi (LTB) model. The LTB solution has recently attracted 
attention, among other things, as a possible explanation of the observed cosmic 
acceleration without resorting to dark energy. Then, we phenomenologically 
constrain both the parameters Ki = —9^/9^ and K2 = —9^ /9^ of the LTB metric 
in the Fermi frame by using different kinds of solar system data. The corrections 
Auj to the standard Newtonian/Einsteinian precessions of the perihelia of the 
inner planets recently estimated with the EPM ephemerides, compared to our 
predictions for them, yield K^ = (4 ± 8) x IQ-^^ s'^, K2 = {3 ± 7) x 10-^3 
s~^. The residuals of the Cassini-based Earth-Saturn range, compared with the 
numerically integrated LTB range signature, allow to obtain Ki p^ K2 ^ 10~^^ 
s~^. The LTB-induced distortions of the orbit of a typical object of the Oort 
cloud with respect to the commonly accepted Newtonian picture, based on the 
observations of the comet showers from that remote region of the solar system, 
point towards Ki ^ K2 ^ 10~^° — 10~^^ s~^. Such figures have to be compared 
with those inferred from cosmological data which are of the order of Ki ^ K2 = 
-4 X 10-36 s-2. 

Subject headings: gravitation — celestial mechanics — astrometry — ephemerides — 
planets and satellites: individual (Mercury, Venus, Earth, Mars) 
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Introduction 



Inhomogeneous cosmological models (JHellabyl |2009[ ) have recently attracted much 
attention because, among other things, some of them may potentially be useful in explaining 
the observed Universe's acceleration without resorting to dark energy. Thus, it is important 
to put them on the test independently of the phenomenon itself for which they have been 
purposely introduced. 

In the framework of the general theory of relativity, the most general form of the line 
element (ds)"^ of a spherically symmetric inhomogeneous space-time in which the source 
in the Einstein's field equations is a perfect fluid is , inu comoving-synchronous spherical 
coordinates {t,r, ^,0} (jPlebahski &: Krasihskil 120061 ) . 



(ds) 



3e:(r,t), 



cdtf 



Mr,t) 



(dr) 



m\r,t)[{deY + sin' eid(i)Y] 
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In order to solve the resulting Einstein field equations, an equation of state has to 
be assumed. The most natural choice consists of setting the pressure equal to zero, 
corresponding to a dust evolution driven by gravitation only. From the equations of motion 
of a perfect fluid it turns out that, for p = 0, the dust moves along timelike geodesies. As a 
consequence, 



€=0 and, foil 9^' ^ 0, 



^ 



l + 2(B{r) 



(2) 



where (£(r) is an arbitrary function such that (B > —1/2 for all r in order to preserve 
the right signature of the metric of eq. ([T]). The resulting space-time line element is 
the so-called Lemaitre-Tolman-Bondi (LTB) model. Indeed, such sp herically symmetric 
in homogeneous du st m odels were fir st discovered by iLemaitrd (119331 ) and further studied 
bv iTolmanI (1934) and iBondil (119471 ): for a detailed discussion of several properties see 
Plebahski fc Krasihskil (120061 ). Such models are among the best known and most useful 
exact solutions of Einstein's equations. Since they allow us to examine non-linear effects 
analytically, or at least in a tractable way, there is an extensive litera ture using them mostly 
to describe cosmological inhomogeneities ( jHellaby fc Krasihskil 120061). b ut also as in other 
theoretical contexts, such as gravitational collapse ( JFranzen et al.ll2010l ) and censorship of 



^Recall that, by definition, the comoving coordinates are those in which the vector field 
m", a = 0, 1, 2, 3 has only the time component, i.e. u"^ (x 5%. If the vector field u" has zero 
rotation, as it happens if the space-time is spherically symmetric and the metric obeys the 
Einstein equations with a perfect fluid source, the comoving coordinates can be chosen so 
that they are synchronous, i.e. in them the metric tensor has no off-diagonal components 
goi, i = 1,2,3. 

^Here and in the following, 9^ denotes the partial derivative of D\ with respect to r. 
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singularities (JDwivedi fc Joshil 119971 ) or quantum gravity flBojowald et al.ll2008[ ). Several 
authors considered LTB models as tools to probe how the cosmic acceleration associated to 
recent observations can be accounted for inhomogeneities, without introducing dark ene rgy 
( IParanjape fc Singhl 120061 : ISarkai] l2008l : lEnqvistl Eoosi : lOarcia-Bellido fc HaugboU^ boOSh . 
LTB models a re also a standard choice to apply Buchert's scalar averaging formalism 
( jBuchertl 120081) . in which the effects of dark energy could be mimicked by "backreaction" 
terms. See ICelerierl ( 120071 ) for a comprehensive review of such aspects. 



Here we are interested in putting constraints in a phenomenological way on some 
parameters of the LTB metric from its local tidal effects on the dynamics of test bodies of 
the solar system. To this aim, it is necessary to set up a standard orthonormal tetrad frame 
parallel transported along the worldline of a fundamental observer representing, in this 
case, the Sun, and explicitly derive the expression of the LTB metric in such a local Fermi 
frame. This has been recently done by iMashhoon et al.l (120071 ). In general, the motion of 
test bodies c an be s t udied i n such Fermi coord i nate systems following i deas and methods 
developed by ISyngd (Il960l ): iManasse k. Misnerl (119631 ): iMashhoonI (119771 ). 



The paper is organized as follows. In Section |2] we analytically work out the LTB effects 
on the orbital motion of a test particle in the local Fermi frame and compare them with 
some results present in literature. In Section [3] we compare our theoretical predictions to 
diff erent types of bodies a nd data of the solar system to phenomenologically put constraints 
on (IMashhoon et al.ll2007| ) 



Ki 






(3) 



Ko 



^ 
dX 



(4) 



entering the LTB metric written in terms of the Fermi coordinates of a fundamental 
observer. Conce r ning their values infered from cosmo logical data, from (62)- (64) by 
Mashhoon et al.l (120071 ) it turns out that, in proximity of the origin of the local Fermi frame, 
i.e. for r — 7- 0, -fCi and K2 are rather similar being of the order ou 



Ki/2 ~ qoHl 



In it (JLarson et al.ll2010l ) 



Hq = 71.0 kms"^ Mpc"^ = 2.3 x 10"^^ s"^ 



is the present value of the Hubble parameter, and ( jXu et al.ll2006l l 

. 1 + 3fiDEWDE 



Qo 



-0.7 



(5) 
(6) 

(7) 



^Ki/2 nieans that both Ki and K2 have approximately the same values. 
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is the deceleration p arameter in which (JLarson et al.ll2010[ ) Qbe = 0.734 ± 0.029 is the dark 
energy density and (JLarson et aLll2010l ) wbe = —^■^'^-oll is the equation of state. Thus, 
the order of magnitude of Ki and K2 inferred from cosmological observations is 



Ki/2 ^ -4 X 10"^^ s-^ 
Section H] is devoted to summarizing our findings and to the conclusions. 



2. Analytical calculation 

In the Fermi coordinates {T,X,Y, Z} associated with a fundamenta l observer 

corre sponding, in this case, to the Sun, the LTB tidal perturbing potentiau (JMashhoon et al. 
2007h 

f/LTB = -I [MX' + Y') + K2Z'] (9) 



arises. It induces a perturbing acceleration ( jMashhoon et al.l 120071 ) 

Altb = -K,R+iK\ - K2)Zk, 
where k is the unit vector along the Z axis. 



(10) 



In order to work out the long-term variations of the Kep lerian orbital el ements, it 
is convenient to adopt the Lagrange's perturbative approach (ICapderoul l2005l ) . valid for 
perturbations arising from a potential function, since, as we will see, it implies just one 
integration. In such a framework, J/ltb has to be evaluated onto the unperturbed Keplerian 
ellipse and averaged over one orbital period, i.e one has to compute 



{Ultb) = ^ / ' Ui^TBdT, 
-Hj Jo 



'ID 



in which [/ltb denotes the perturbing potential computed in terms of the parameters of 
the Keplerian orbit. To this aim, useful relations in terms of the eccentric anomaly E are 



^Here and in the following application of the Lagrange's equations for the variation of 
the Keplerian orbital elements U is define d according to the conv ention A = Vf/. That is 



why our ?7ltb of eq. (|9]) is Vn of (A.l) by iMashhoon et al.l (120071 ) with a minus sign. 
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f lCapderoull2005f ) 



R 



a(l — ecos-E), 



(12) 



cos E — e 

cos/ = -, 

1 — e cos h 



(13) 



sin/ 



dT 



a/1 — e'^sinE 
1 — e cos E 

(1 — ecosi?) 



rfE, 



n 



(14) 



(15) 



where a and e are the semimajor axis and the eccentricity, respectively, of the unperturbed 
Keplerian elhpse, / is the true anomaly reckoning the position of the test particle along it, 
and n = 2tx / P\, = ^JCMja^ is the unperturbed Keplerian mean motion. The Cartesian 
coordinates of the m oving test partic le along the unperturbed Keplerian ellipse, explicitly 
entering eq. (|9]), are (ICapderoull2005l ) 



X = R (cos r2 cos u — cos / sin f2 sin xi) 



(16) 



Y = R (sin Q cos u + cos / cos Q sin u) 



(17) 



i? sin J sin M, 



where Q is the longitude of the ascending node, / is the inclination of the orbit to the 
reference {^1^} plane, u = u + f is the argument of latitude, in which u is the argument of 
pericentre. Thus, the perturbing potential, averaged over one orbital revolution, turns out 
to be 

(Uutb) = K{V{a, e, /, u) + i^2^(a, e, /, w), (19) 



with 



and 



V 



3 + cos 2H (9 + 3 cos 2/ + 10 sin^ / cos 2uj) 



£ 



a? sin I 

4 



1- — (-3 + 5 cos 2a;) 



(20) 



(21) 



Now, eq. f ITI?]) . with eq. f l2U]) and eq. f l?I]) . can be plunged into the right-hand-sides of the 
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Lagrange's equations which are f ICapderoul |2005| ) 

da 
dT 



na \ dM 



de 
dT 

dl 

dT 

dVL 

It 

dcu 
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1 1 



na^ 



1 dU dU \ 



1 



na^ a/I — e^ sin / 
1 



_dU_ 






na^y/l — e^ sin/ \dl J 

1 n-e^dU cosIdU 



na- 



■■VT^ 



n + 



1 



na 



dU 
^ \ da 



de sin/ dl 

l-e^dU 



de 



(22) 
(23) 
(24) 
(25) 
(26) 
(27) 



where M. is the mean anomaly. In this way, the long-term variations of the Keplerian orbital 
elements can be straightforwardly obtained from simple derivatives. From eq. (1221) and eq. 
( TT9l) . with eq. ( l20l) and eq. ( 12T|) . it immediately turns out that a remains unchanged, on 
average. Instead, the other Keplerian orbital elements undergo long-term changes. They 
are 



/ = 



5e\/l - e'^AK . , r • <. 

sm ismzoj, 

An 

= sin 2/ sin 2a;, 



n) = 



8nVl - e^ 
A/i cos / 



1 + — (3-5cos2a;) 



7W) = 



2nVl - e2 

(-11 + 6e2)/i:i + (-1 + 6e'^)K2 + 5AK [(-1 + 2e2)cos2a; - 2cos2/sin2 w] 

8nVl - e2 

(7 + 3e2) rs/sTi + /sTa + A-fsT cos 2/ + 10(1 + e2) Aii' sin^ I cos 2w] 

"+ 8^ ' 



(28) 
(29) 
(30) 
(31) 
(32) 



where we have defined AK = Ki — K2. While the eccentricity e and the inclination / 
experience only long-term, harmonic variations with frequency 2u, the longitude of the 
ascending node Q, the argument of pericenter uj and the mean anomaly Ai undergo secular 
precessions as well. 

Much more cumbersome calculations show that the Gausqfl perturbative approach 
( ICapderoul 120051 ) . based on the projections of the perturbing acceleration Altb of eq. flTOj) 



^Contrary to the Lagrange scheme, it is not limited to accelerations derivable from 



onto the radial R, transverse H and normal Y directions of a frame co-moving with the 
test particle and implying six integrals, yields the same results as eq. ( l28l) -eq. ( !32l) . and 
also (a) = 0. 

It turns out from eq. (l30!) -eq. (!3T|) that the precession of the longitude of pericentre 
w = Q + u can be cast into the useful form 



(ttr) = KiT{I, e, uj) + K2G{I, e, oo), 



with 
and 



■^ + cos/ 



I cos 2/ 



I sin^ / cos 2w + e^ [|(l + cos/) + 5sin2 (|) cos 2a;] 



sm' 



(I) 



2nVT 



2n^ 



[3 — 5 cos 2a; + 10 cos / sin^ a; + e^ (5 cos 2a; — 3 



(33) 



(34) 



(35) 



From eq. ( J34l) and eq. ( J35l) it can be noted that both secular and long-period, harmonic 
terms are present. For orbits exhibiting moderate eccentricities and inclinations to the 
reference plane, like those of the planets of the solar system, the harmonic terms play a 
minor role because they appear multiplied by sin^/, sin^(J/2), e^sin^(//2). Anyway, since 
we do not know a-priori the magnitudes of K^ and K2^ we will retain all the terms of eq. 
f p^ and eq. fl35|l in our calculations. 



Concerning the results by iMashhoon et al.l (120071 ) . they used the Gauss equat i ons fo r 
the variation of the elements and the Delauney orbital elements. iMashhoon et al.l (120071 ). 
interested in slightly eccentric orbits, obtained non-zero secular variations of the semimajor 
axis a of order 0[e^\ and of the eccentricity e of order 0[e\ although not explicitly shown; 
actually, eq. (|28|l is of order 0(e), w hile we obtain (a) = to all orders in e since it is an 
exact result. iMashhoon et al.l (120071 1 displayed the secular precessions of the argument of 
pericentre a;, to order 0(e), and of th e longitude of th e ascen ding node fi, to order C(e^). 
While eq. (|30|) agrees with (A. 25) by IMashhoon et al.l (120071 ). eq. (|3T]) is in disagreement 
with the sum of ( A. 24) for a; and the additional pericentre precession due to A = —KiR by 



Mashhoon et al.l ( 2007). even in the limi t e — ?• 0. The variation of the mean anomaly was 



Mashhoon et al.l ( 120071 ). 



not computed by 

It is also interesting to work out the variations after one orbital revolution of the radiaJ, 
transverse and normal components of the perturbation of the radius vector R. By usingj 



potential function, being valid for any kind of perturbations. 

^Here the shifts in the Keplerian orbital elements come from indefinite integration(s), so 
that they are functions of the fast angular variable used which, in our case, is E. 
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f lCasottolll993h 



AR 



- Aa - a cos /Ae + ae(l - e^y^/^ sin fAM 



A^ = a sill / 



a(l-e2) 

AT = 7' (sinuA/ — cosusin/Ari) , 
it is possible to obtaiiJll 

bnaey/l - e'^AK sin^ / sin 2uj 



Ae + r (Aw + cos /AO) H \/\ - e^ AM 



AR = 



AE 



AT = 



2n2 



1 + e 7ra { (2 + 3e) [3A'i + /4:2 + AT^f cos 21] + WeAK sin^ / cos 2w} 



1-e 



2n2 



i(-l + e)Vl - e2Ai4'sin2/cosi: 
2^^2 



(36) 

(37) 
(38) 

(39) 
(40) 
(41) 



Concerning the perturbation of the vplop ity vpptnr ^ its radial, transverse and normal 
components can be expressed, in general, as (ICasottolll993l ) 



AVr 

AVs 
AW 



n sin / f e 



, , -Aa-\ Ae 

VT^l? \2 r 



na^\pl 



(Aw + cos/Af}) 5-AA^, 



ivT 



2r 



-Aa 



na(e + cos /) nae sin / 



(1 - e2)3/2 



Ae 



yr^^ 



(Aw + cos/Af^), 



x/T^^ 



[(cos u + e cos w) A/ + (sin u + e sin w) sin / Afi] . 



From them it is possible to obtai 

7ra { [-4 - 3e(l + 2e - e2)] {iKi + K2 + AK cos 21] + 10e(-l - 2e + e^)AK sin^ I cos 2a;} 



AVp 



AVh 



AW 



4(l-e)2n 



BiraeAK sin / sin 2a; 
2(1 -e)n ' 

7ra(l + Ae^)AK cos / sin 2a; 
2(1 -e)n ■ 



(42) 

(43) 
(44) 

{,45) 
(46) 
(47) 



While the transverse and normal components experience only long-term harmonic 
variations, the radial one exhibits a secular variation as well. 



^Here A^ = Ae(27r) - Ae(0), ^ = R,E, T. 

^Also in this case, it is to be intended AVj = AVj{2tt) - AV,(0), j = R, S, T. 
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3. Confrontation with the latest observational determinations 
3.1. The precessions of the perihelia of the inner planets of the solar system 



Recently, iPitjeval ( l2010l ) has analyzed more than 550000 planetary observations of 
several kinds covering the time interval 1913 — 2008. She used the dynamical force models 
of the EPM2008 ephemerides by estimating about 260 parameters. Among them, she also 
determined corrections Arij to the standard Newtonian/Einsteinain perihelion precessions of 
all the planets of the solar system including Pluto as well; such corrections, by construction, 
account for any unmodelled/mismodelled dynamical effects, so that they can be used, in 
principle, to put constraints on Ki and K2. To this aim, we will use the inner planets, 
whose estimated perihelion corrections are listed in Table [H because they are more accurate. 
In Table [2] we compute the values of the coefficients J-" and Q, in s, for the inner planets 
of the solar system. By equating the predicted precession of eq. (155]) . with the figures of 
Table [21 to the estimated corrections Aw listed in Table [H it is possible to write down six 
inhomogeneous linear systems of two equations in the two unknowns Ki and K2. It turns 
out that the Earth and Mars yield the tightest constraints on Ki and K2: they are 

Ki = (-1.2 ± 1.4) X 10"^^ s"^ (48) 

K2 = (1.4±1.5) X 10"^^ s"l (49) 

The quoted uncertainties have been obtained by linearly propagating the errors in Aw of 
Table [H 



In Table |3] we quote the values of the corrections Aw estimated by iPitjeval ( 120051 ) with 
older versions of the EPM ephemerides and less extended data sets; for Venus, Earth and 
Mars they are, for some reasons, more accurate by about one order of magnitude than the 
more recent results of Table [U Using the values of the corrections Aw of Table [3] allows to 
obtain more stringent constraints for Ki and K2. It turns out that, in this case, Venus and 
the Earth yield the tightest ones which are 

Ki = (4 ± 8) X lO-^S-^ (50) 

K2 = (3 ± 7) X 10"^^ s~^ (51) 



Table 1: Estimated corrections Aw, in mas cty ^ (1 mas cty -"^ = 1.5 x 10 ^^ s ^), to the 
standard perihelion precessions with the EPM2008 ephemeri des. The quote d errors are not 



the formal, statistical ones but are realistic. From Table 8 of IPitjeval ( 120101 ) 



Mercury Venus Earth Mars 



Arij (mas cty^i) -4 ± 5 24 ± 33 6±7 -7 ±7 
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In both cases, it turns out that Ki,K2 and AK are statistically compatible with zero; 
anyway, note that K2 may be up to three orders magnitude larger than Ki. 

It may be interesting to use such values for Ki and K2, derived from the inner planets, 
with Saturn in order to see if they yield a LTB perihelion precession for the ringed planet 
compatible with its anomalous perihelion precession 



Adj^^f ^^ = -0.006 ± 0.002 arcsec cty-\ (52) 

Atij^^f ") = -0.010 ± 0.015 arcsec ctj-\ (53) 



Avu^^^^> = -0.010 ± 0.008 arcsec cty"^ (54) 

recently estimated by processing also some years of normal points from the Ca s sini 
space craft with both thqj EPM and the INPOP ephemerides ( Fienga et al.ll2010l : IPitjeva 



2OIOI ): see Section [3^ for an analysis involving a different quantity related to such 



interplanetary ranging data. To this aim, we have, first, to compute 

J'sat = -2.21546 X 10^ s, (55) 



Then, it turns out that 



^sat = -22134.71s. (56) 



(5(Ali7Sat) + |J^Sat|5i^l + |^Sat|5A'2 



< 1 (57) 



for all the values of AriJsat (eq. flS^ -eq. (JSID) and the previously obtained figures for Ki 
and K2 (eq. (gHD-eq. (gg]) and eq. ([50D-eq. (|5T|)). 
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The value of eq. fl52|) has been communi cated by E.V. Pitjeva to the author in November 



2008, and it has explicitly been reported by iFienga et al.l (120 lOf ) in their Table 4. 



Table 2: Numerical values, in s, of the coefficients J-" and Q for the inner planets according 
to eq. flM|) and eq. f l35|) . The semimajor axis a, entering the Keplerian mean motion n, 
the eccentricity e, the inclination / and the a rgument of perihelion u have been retrieved 



from Table A. 2 of iMurray and DermottI ( 119991 ). Reference frame: ICRF/J2000. Coordinate 



system: Ecliptic and Mean Equinox of Reference Epoch. 



Mercury Venus Earth Mars 

F (s) -1.76948 X 10« -4.61929 x 10^ -7.53291 x 10^ -1.40874 x 10^ 
g (s) -6238.21 -15397.9 -6.99848 x 10^^ -20198.2 
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3.2. The Earth-Saturn range 



Actually, it is possible to infer lower upper bounds by using different targets. After the 
Cassini spacecraft started its "grand tour" of the Saturnian system, it has been possible to 
drastically increase the accuracy of the orbit de termination of the ri nged planet through 



direct ranging to Cassini itself. Figure B-20 of iFolkner et al.l ( l2008l ) shows the range 
residuals of Saturn from 2004 to 2006 constructed with the DE421 ephemerides from 
Cassini normal points; processing of extended data reco rds of Cassini is currently ongoing. 



The range residuals of Figure B-20 (IFolkner et al.ll2008l ) are accurate at 10 m level. Also a 



pair of range residuals from close encounters with Voyager 1 (1980) and Voyager 2 (1982) 
are shown: they are almost one order of magnitude less accurate. Thus, we will, now, 
look at the perturbation induced by Ki and K2 on the Earth-Saturn range p and we 
will compare it to the DE-421-based residuals, which are a somewhat more direct kind of 
observable with respect to the corrections to the perihelion precession^lj. Figure [1] shows 
the numerically computed LTB Earth-Saturn range perturbation for the largest values of 
eq. flSU]) -eq. flHTl) over a time interval of 2 yr, comparable to that of the available Cassini 
residuals. The amplitude of the signal is 1750 m, which is two orders of magnitude larger 
than the Cassini residuals. Even by assuming that part of the putative LTB signal, not 
explicitly modeled in the dynamical force models used to construct the residuals, has been 
absorbed and removed in the estimation of the initial conditions, it seems reasonable to 
conclude that a non-negligible part of a so huge anomalous effect would have been left in 
the residuals and, thus, it would have not escaped from detection. It can be shown that 
by posing Ki = K2 ^ 10~^^ s~^ yields a LTB signal as large as 100 m over 2 yr, while for 
Ki = K2 ~ 10"^^ s~^ the amplitude of the LTB range is just 10 m over 2 yr. Thus, we 
conclude that Saturn poses 10~^^ s~^ as upper bounds for both Ki and K2. It is likely that 
when the analysis of more extended Cassini data sets will be completed, more stringent 
limits on Ki and K2 will be placed. 



^°The longitude of the perihelion, as all the other Keplerian orbital elements, are not 
directly observable quantities. 



Table 3: Corrections Avj, in milliarcsec cty~^ (1 mas cty^^ = 1.5 x 10^^^ s'^), to the standard 
perihelion precessions estimated by E.V. Pitjeva with the EPM2005 (Mercury, Earth, Mars) 
and EPM2006 (Venus) ephemerides . The quoted errors are not the formal, statistical ones 
but are reali s tic. F rom Table 3 of IPitjeval f 20051 ) (Mercury, Earth, Mars) and Table 4 of 



Fienga et al.l f l2010f ) (Venus 



Mercury Venus Earth Mars 



Aro (mas cty^i) -3.6 ± 5.0 -0.4 ± 0.5 -0.2 ± 0.4 0.1 ±0.5 
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Earth-Saturn range: SSB numerical calculation 
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Fig. 1. — Difference A|p| in the numerically integrated (Earth-Moon-Barycenter)-Saturn 
ranges with and without the nominal perturbation due to the LTB cosmological tidal effects 
for Ki = 12 X 10~^^ s~^ and K2 = 7 x 10~^^ s~^ over AT = 2 yr. The same initial conditions 
(J2000) have been used for both the integrations. The state vectors at the reference epoch 
have been retrieved from the NASA JPL Horizons system. The integrations have been 
performed in the ICRF/J2000.0 reference frame centered at the Solar System Barycenter 
(SSB). 



3.3. The Oort cloud 



The Oort cloud (jOortlll950l ). populated by a huge number of small bodies moving 
along very eccentric orbits highly inclined to the ecliptic, is supposed to exist in the remote 
peripheries of the solar system. Indeed, it is a spheroid with a semimajor axis of about 100 
kau e longated towards the Galactic center and a semiminor axis of about 80 kau (jWeissman 
19961 ). The interaction of a nearby passing star with the Port comet cloud can give rise to 
comet showers reaching the region of the major planets (lHillslll98ll : IWeissmanI Il996[ ) . The 
formation time scale of such a comet shower is about 1 Myr. For example, according to 
BobvlevI fl2010f ). the star HIP 89 825 (GL 710) has a probability of 0.86 of penetrating the 
Oort cloud in the next 1.45 ± 0.06 Myr. Moreover, its has a smaller, non-zero probability 
(1 X 10~^) of falling into the region r < 1 kau, thus potentially influencing the dynamics of 
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the Edgeworth-Kuiper belt objects as well. 

It is interesting to see the extent to which the LTB cosmological tidal effects may 
alter the standard picture of the Oort cloud. To this aim, we will numerically integrate 
over a (Keplerian) orbital period the paths of a typical Oort comelo {Xq = 40 kau, 
Yq = 30 kau, Zq = 5 kau, Xq = —23 kau Myr~^, Yq = —15 kau Myr~^, Zq = —15 kau 
Myr~^ corresponding to a = 66.7 kau, e = 0.92, / = 81.6 deg), sharing the same initial 
conditions, with and without the LTB tidal acceleration of eq. f llOp . Although not in a 
strict quantitative way as in the previous cases examined so far, such an approach may 
indicatively provide us with tighter constraints on Ki and A'2 by excluding those values 
yielding too bizarre or implausible scenarios for the orbits of the Oort comets. 

Actually, Figure |2l depicting the sections in the coordinate planes of the purely 
Newtonian and LTB orbits for Ki p^ K2 ~ 10~^^ s~^, and Figure [31 showing the heliocentric 
distances and velocities of the same paths, tell us that the values for the LTB parameters 
Ki and K2 allowed by the Earth-Saturn range are likely still too large. Indeed, the 
resulting numerically integrated paths are radically different from the Newtonian ones. 
It must also be noted that the LTB acceleration (3.0 x 10"^^ m s"^) is larger than the 
Newtonian one (2 x 10~^^ m s~^), so that the perturbative scheme would not be applicable: 
the numerical integration of the equations of motion including both the LTB and the 
Newtonian accelerations is, thus, strictly required. The consequences on the dynamics of 
the Oort cloud and on the comet showers may be not negligible, although, of course, not 
easily predictable. Suffices it to note from Figure |3] the globally smaller spatial extension 
of the LTB orbit and its higher velocity which might reduce the impact of nearby passing 
staro. Further numerical analyses show that reasonable orbits, i.e. not too far from the 
Newtonian picture, occur for ^1/2 < 10"^° - 10"^^ s"^, when ALTBMNcwton ~ 10"^ - 10"^. 

It must be stressed that it turns out that, actually, the LTB trajectories exhibit a 
strong dependence on the initial conditions. Thus, a complete statistical analysis involving 
vast ensembles of different initial conditions would be required, but it is outside the scope 
of the present work. 



^^It can be shown that, in the Newtonian case, such initial conditions yield X^ax = 90 
kau, Fmax = 60 kau, ^max = 80 kau over on e Keplerian orbital period, so that it is well within 



the domain indicated by IWeissmanI (119961 ). 

^^Clearly, also the motion of a perturbing star should be worked out in the framework of 
the LTB dynamics, but this is beyond the scope of the present paper. 
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Fig. 2. — Sections in the coordinate planes of the numerically integrated orbits of an Oort 
comet with a = 66.7 kau, e = 0.92, / = 81.6 deg. Dashed blue line: Newton. Dash-dotted 
red line: Newton-(-LTB with Ki = A x 10"^^ s"^ A'2 = 3 x 10"^^ s~^ The initial conditions 
are Xq = 40 kau, Yq = 30 kau, Zq = 5 kau, Xq = —23 kau Myr"-'^, Yq = —15 kau Myr~^, Zq = 
— 15 kau Myr~^ (1 kau Myr~^ = 0.0047 km s~^). The time span of the integration is 
-Ph<T < 0. 
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Fig. 3. — Heliocentric distances, in kau, and velocities, in km s~^, for the numerically 
integrated orbits of an Oort comet with a = 66.7 kau, e = 0.92, I = 81.6 deg. Dashed 
blue line: Newton. Dash-dotted red line: Newton -l-LTB with Ki = A x 10^^^ s^^, K2 = 
3 X 10"^^ s~^. The initial conditions are Xq = 40 kau, Yq = 30 kau, Zq = 5 kau, Xq = 
-23 kau Myr-\ % = -15 kau Myr-\ Zq = -15 kau Myr'^ (1 kau Myr'^ = 0.0047 km s^^). 
The time span of the integration is — Pb < 7" < 0. 
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4. Summary and conclusions 

We, first, analytically worked out the long-term, i.e. averaged over one orbital 
revolution, perturbations on the orbit of a test particle moving in a local Fermi frame 
induced by the cosmological tidal effects of the inhomogeneous Lemaitre-Tolman-Bondi 
model which has recently attracted attention as possible explanation of the observed cosmic 
acceleration without resorting to dark energy. In particular, we computed the variations of 
the semimajor axis a, the eccentricity e, the inclination J, the longitude of the ascending 
node Q, the argumet of pericenter uj and the mean anomaly Ai by means of the Lagrange's 
variational equations with the eccentric anomaly E as fast variable. We found that a 
does not experience any change, on average, while the other Keplerian orbital elements 
undergo long-term variations including both secular and harmonic terms with frequency 
2u] e and / exhibit only sinusoidal changes, while I,Q,Ai secularly precess as well. We 
repeated the calculation also with the Gauss perturbative approach by finding the same 
results. We also explicitly computed the changes over one orbital revolution of the position 
and velocity vectors along the radial R, transverse S and normal directions T of a frame 
co-moving with the test particle. While the radial and normal components AR and AT of 
the perturbation of the position vector experience only harmonic variations, the transverse 
one AS shows secular changes as well. Concerning the velocity, its radial perturbation 
AVr undergoes both secular and sinusoidal modifications, while its transverse and normal 
components AV= and AVf exhibit harmonic signatures only. In obtaining our results we 
made no approximations on e and I by retaining all terms. 

Then, we phenomenologically put constraints on both the parameters Ki and K2 
of the LTB metric in the local Fermi frame by looking at various astronomical bodies 
and data of the solar system. By comparing our analytical prediction for the rate of 
the longitude of pericenter vu to different sets of the corrections Aw to the standard 
Newtonian/Einsteinian precessions of the perihelia of the inner planets recently estimated 
with the EPM ephemerides we found that the tightest constraints are Ki = {A ± 8) x 10^^^ 
s~^, i^2 = (3 ± 7) X 10~^^ s"^. The confrontation of the residuals of the Earth-Saturn range, 
obtained by processing some years of radiotechnical data from the Cassini spacecraft as 
well, with the numerically computed LTB-induced Earth-Saturn range signal allowed to 
set Ki ^ K2 ^ 10~^^ s~^. By looking at the LTB-type distortions of the orbit of a typical 
object of the Oort cloud with respect to the commonly accepted Newtonian picture, based 
on the observations of the comet showers from that remote region of the solar system, 
pointed towards even lower bounds, i.e. Ki ^ K2 ^ 10"^" — 10~^^ s~^. According to 
cosmological data, Ki ^ K2 = —4 x 10~^® s~^. 
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